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Abstract. We propose an interpolation formula for the distribution of the reflection 
coefficient in the presence of time reversal symmetry for chaotic cavities with 
absorption. This is done assuming a similar functional form as that when time reversal 
invariance is absent. The interpolation formula reduces to the analytical expressions 
for the strong and weak absorption limits. Our proposal is compared to the quite 
complicated exact result existing in the literature. 
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1. Introduction 

In recent years there has been a great interest in the study of absorption effects on 
transport properties of classically chaotic cavities [DElElllElinilZllHliniliniiniliaiini 

HU 113 Cni im CH] (for a review see Ref. jUI). This is due to the fact that for experiments 
in microwave cavities |2IIl E] , elastic resonators j22] and elastic media [SB] absorption 
always is present. Although the external parameters are particularly easy to control, 
absorption, due to power loss in the volume of the device used in the experiments, is an 
ingredient that has to be taken into account in the verification of the random matrix 
theory (RMT) predictions. 

In a microwave experiment of a ballistic chaotic cavity connected to a waveguide 
supporting one propagating mode, Doron et al [T] studied the effect of absorption on 
the 1x1 sub-unitary scattering matrix S*, parametrized as 



where R is the reflection coefficient and 6 is twice the phase shift. The experimental 
results were explained by Lewenkopf et al. j2| by simulating the absorption in terms 
of Np equivalent "parasitic channels", not directly accessible to experiment, each one 
having an imperfect coupling to the cavity described by the transmission coefficient Tp. 
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A simple model to describe chaotic scattering including absorption was proposed 
by Kogan et al. 0. It describes the system through a sub-unitary scattering matrix 
S, whose statistical distribution satisfies a maximum information-entropy criterion. 
Unfortunately the model turns out to be valid only in the strong-absorption limit and 
for i? <^ 1. For the 1x1 S-matrix of Eq. ([Q), it was shown that in this limit 9 is 
uniformly distributed between and 27r, while R satisfies Rayleigh's distribution 

P0{R) = ae-"^; < 1, and a > 1, (2) 

where /? denotes the universality class of S introduced by Dyson ^j: /? = 1 when time 
reversal invariance (TRI) is present (also called the orthogonal case) , /? = 2 when TRI 
is broken {unitary case) and /? = 4 corresponds to the symplectic case. Here, a = 7/3/2 
and 7 = 27r/raA, is the ratio of the mean dwell time inside the cavity (27r/A), where A 
is the mean level spacing, and is the absorption time. This ratio is a measure of the 
absorption strength. Eq. Q is valid for 7 3> 1 and for i? <C 1 as we shall see below. 

The weak absorption limit (7 ^ 1) of PpiR) was calculated by Beenakker and 
Brouwer [Hj, by relating R to the time-delay in a chaotic cavity which is distributed 
according to the Laguerre ensemble. The distribution of the refiexion coefficient in this 
case is 

In the whole range of 7, PpiR) was explicitly obtained for /? = 2 j3]: 

-7/(l-R) 

P2iR) = - 1) + (1 + 7 - e^)(l - R)] , (4) 

and for /5 = 4 more recently [12]. Eq. (jlj reduces to Eq. Q for small absorption (7 <^ 1) 
while for strong absorption it becomes 

Notice that P2{R) approaches zero for R close to one. Then the Rayleigh distribution, 
Eq. (j21), is only reproduced in the range of few standard deviations i.e., for R ~ 7"^. 
This can be seen in Fig. ^a) where we compare the distribution P2{R) given by Eqs. 
and (jSj) with the exact result given by Eq. (0} for 7 = 20. As can be seen the result 
obtained from the time-delay agrees with the exact result but the Rayleigh distribution 
is only valid for i? <C 1. 

Since the majority of the experiments with absorption are performed with TRI 
(/? = 1) it is very important to have the result in this case. Due to the lack of an exact 
expression at that time. Savin and Sommers |H] proposed an approximate distribution 
Pj3=i{R) by replacing 7 by 7/3/2 in Eq. (j3)). However, this is valid for the intermediate 
and strong absorption limits only. Another formula was proposed in Ref. as an 
interpolation between the strong and weak absorption limits assuming a quite similar 
expression as the (3 = 2 case (see also Ref. [TB])- More recently [17], a formula for the 
integrated probability distribution of x = (1 + -R)/(1 — i?), W{x) = PQ^~^\x)dx, was 
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obtained. The distribution P^=i(i?) = ji^^Pq (t^) then yields a quite comphcated 
formula. 

Due to the importance to have an "easy to use" formula for the time reversal case, 
our purpose is to propose a better interpolation formula for P/3{R) when (3 = 1. In the 
next section we do this following the same procedure as in Ref. ^H]- We verify later 
on that our proposal reaches both limits of strong and weak absorption. In Sec. El we 
compare our interpolation formula with the exact result of Ref. ^7] . A brief conclusion 
follows. 



2. An interpolation formula for /3 = 1 

From Eqs. Q and we note that 7 enters in Pp{R) always in the combination 7/3/2. 
We take this into account and combine it with the general form of ^(-R) and the 
interpolation proposed in Ref. For /? = 1 we then propose the following formula 

for the i?-distribution 



a 



i/2(e"_i) + (i + a-e")2Fi 



1 1 „\ 1 - 



2'2' ''7 2 



where a = 7/2, 2P1 is a hyper-geometric function ^25j, and Ci{a) is a normalization 
constant 

^'^"^ " (e" - l)r(3/2, a) + aV2(i + a- e")/(a)/2 



where 



/(«)=/ Si^fJII,!;!--] (8) 



la ' '\2'2 X 

and r(a,x) is the incomplete F-function 

00 



F(a,x) = / e-'f-'dt. (9) 

In the next sections, we verify that in the limits of strong and weak absorption we 
reproduce Eqs. and Q. 



3. Strong absorption limit 

In the strong absorption limit, a 00, F(3/2,q;) a^^^e"", and /(a) a;~^/^e~". 
Then, 

lim CUa) = r-nT--7 rr- - (10) 

^ ^ (e°-l)ai/2 + (l + a-e")/2 ^ ^ 

Therefore, the /^-distribution in this limit reduces to 

Q,g-QR/(l-ii) 

^l(^) ^ (1 _ ^)5/2 " » 1' (11) 

which is the equivalent of Eq. but now for (3 = 1. As for the (3 = 2 symmetry, it is 
consistent with the fact that Pi{R) approaches zero as R tends to one. It reproduces 
also Eq. (0) in the range of a few standard deviations {R ~ 7^^ ^ 1), as can be seen in 
Fig.Ottb). 
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Figure 1. Distribution of the reflection coefficient for absorption strength 7 = 20, 
for (a) /3 = 2 (unitary case) and (b) /3 = 1 (orthogonal case). In (a) the continuous 
line is the exact result Eq. while in (b) it corresponds to the interpolation formula, 
Eq. ||HJl. The triangles in (a) are the results given by Eq. Q for /3 = 2 and in (b) they 
correspond to Eq. Hll|) . The dashed line is the Rayleigh distribution Eq. , valid only 
for R ~ 7~-^ and 7 3> 1. 



4. Weak absorption limit 

For weak absorption a —>■ 0, the incomplete F-function in Ci{a) reduces to a F-function 
F(x) [see Eq. Q]. Then, Pi{R) can be written as 



a + «72 + ■ ■ ■)r(3/2) - («V2/2 + . . ■)/(0)/2 



- (aV2 + a^Q + ■ ■ •)2i^i(l/2, 1/2, 1; R){1 - R)/2]. (12) 
By keeping the dominant term for small a, Eq. Q is reproduced. 

5. Comparison with the exact result 

In Fig.|21we compare our interpolation formula, Eq. with the exact result of Ref. pTj . 
For the same parameters used in that reference we observe an excellent agreement. In 
Fig. 121 we plot the difference between the exact and the interpolation formulas for the 
same values of 7 as in Fig. |21 The error of the interpolation formula is less than 4%. 



6. Conclusions 



We have introduced a new interpolation formula for the reflection coefficient distribution 
P/3{R) in the presence of time reversal symmetry for chaotic cavities with absorption. 
The interpolation formula reduces to the analytical expressions for the strong and weak 
absorption limits. Our proposal is to produce an "easy to use" formula that differs 
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Figure 2. Distribution of the reflection coefficient in the presence of time- reversal 
symmetry for absorption strength 7=1,2, 5, and 7. The continuous lines correspond 
to the distribution given by Eq. (jnj. For comparison we include the exact results of 
Ref. [ni (dashed lines). 



by a few percent from the exact, but quite complicated, result of Ref. ^7j. We can 
summarize the results for both symmetries (/? = 1, 2) as follows 



-a/{l-R) 



l) + (l + a 



[1 - Rf+PI^ 

where C/3(q;) is a normalization constant that depends on a = 7/3/2. This interpolation 
formula is exact for (3 = 2 and yields the correct limits of strong and weak absorption. 



(13) 
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